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Abstract. We study the properties of chaotic motions in the intra cluster medium using a set of 20 galaxy clusters simulated with 
large dynamical range, using the Adaptive Mesh Refinement code ENZO (e.g. Norman et al.2007). The adopted setup allows 
us to study the spectral and spatial properties of turbulent motions in galaxy clusters with unprecedented detail, achieving an 
maximum available Reynolds number of the order of R e ~ 500 - 1000 for the largest eddies. The correlations between the 
energy of these motions in the Intra Cluster Medium and the dynamical state of the host systems are studied, and the statistical 
properties of turbulent motions and their evolution with time support that major merger events are responsible for the injection 
of the bulk of turbulent kinetic energy inside cluster. Turbulence is found to account for a ~ 20 - 30 per cent of the thermal 
energy in merging clusters, while it accounts for a ~ 5 per cent in relaxed clusters. A comparison of the energies of turbulence 
and motions in our simulated clusters with present upper-limits in real nearby clusters, recently derived with XMM-Newton 
(Sanders et al.2010), is provided. When the same spatial scales of turbulent motions are compared, the data from simulations 
result well within the range presently allowed by observations. Finally, we comment on the possibility that turbulence may 
accelerate relativistic particles leading to the formation of giant radio halos in turbulent (merging) clusters. Based on our 
simulations we confirm previous semi-analytical studies that suggest that the fraction of turbulent clusters is consistent with 
that of clusters hosting radio halos. 
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1. Introduction 

It is generally believed that turbulent motions are generated 
during the hierarchical assembly of matter in galaxy clusters 
that occurs through mergers and accretion of satellites. 

Mergers between galaxy clusters are very energetic events 
where a fraction of the kinetic energy of Dark Matter (DM) 
halos is transferred into the thermal and kinetic energy. The 
gas from infalling halos is stripped within cluster cores due to 
the action of ram pressure and several fluid instabilities. This 
drives turbulent motions, which may eventually transfer en- 
ergy from large to smaller scales (e.g. Sarazin 2002; Cassano 
& Brunetti 2005; Subramanian et al. 2006; Brunetti & Lazarian 
2007). Also the sloshing motions of DM cuspy cores that may 
occur in the cores of clusters at later stages of merging events 
can drive turbulent motions (e.g. Markevitch & Vikhlinin 2007; 
Roediger et al.2010). Even without considering the process 
of hierarchical formation, galaxy clusters host several poten- 
tial sources of turbulent motions in the Intra Cluster Medium 
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(ICM). The activity of AGN that are found in the central cluster 
galaxies may excite turbulence in the cores of the hosting clus- 
ters (e.g. Fabian et al. 2003; Rebusco et al.2005; Scannapieco 
& Briiggen 2008), as well as galaxies that cross the cluster at- 
mosphere (e.g. Deiss et al 1996; Subramanian et al. 2006). 

From the theoretical viewpoint, turbulence may have a 
deep impact on the physics of the ICM (e.g. Narayan & 
Medvedev 2001; Shekochihin et al 2005, 2010; Lazarian 2006; 
Ruszkowsky & Peng 2010) and on the properties of non- 
thermal components in galaxy clusters (e.g. Subramanian et al. 
2006; Brunetti & Lazarian 2007; Brunetti & Lazarian 2010). 
The presence of turbulent gas motions in the ICM is suggested 
by measurements of the Faraday Rotation of the polarization 
angle of the synchrotron emission from cluster radio galaxies. 
These studies show that the magnetic field in the ICM is tan- 
gled on a broad range of spatial scales (e.g. Murgia et al. 2004; 
Vogt & Ensslin 2005; Bonafede et al. 2010; Govoni et al. 2010; 
Vacca et al. 2010) suggesting the presence of super-Alfvenic 
motions in the medium. Also independent attempts from X-ray 
observations of a number of nearby clusters, based on pseudo- 
pressure maps of cluster cores and on the lack of evidence for 
resonant scattering effects in the X-ray spectra, provided hints 
of turbulence in the ICM (e.g. Schuecker et al. 2004; Henry et 
al. 2004; Churazov et al. 2004; Ota et al. 2007). Important con- 
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Fig. 1. 2-dimensional slice showing the gas temperature for the innermost region of galaxy cluster El, during its main merger 
event (z = 0.6). The side of the slice is 8.8 Mpc h and the depth along the line of sight is 25 kpc ft -1 . 



straints on the fraction of the turbulent and thermal energy in 
the cores of clusters are based on the analysis of the broadening 
of the lines in the emitted X-ray spectra of cool core clusters 
(Churazov et al.2008; Sanders et al. 2010a,b). In several cases 
these studies derived interesting upper limits to the ratio of tur- 
bulent and thermal energy, of the order of < 20 per cent, even 
if the most stringent limits refer to clusters with very compact 
cores and are thus sensitive to (the energy of) motions on a few 
tens kpc scale. 



Nowadays numerical simulations provide a unique way to 
study the generation and evolution of velocity fields in the ICM 
across a wide range of scales. Early Eulerian numerical sim- 
ulations of merging clusters (e.g., Roettiger, Stone & Burns 
1999; Norman & Bryan 1999; Ricker & Sarazin 2001) ob- 
tained the first reliable representation of the way in which tur- 
bulence is injected into the ICM by merger events. More re- 
cent works, performing high-resolution Lagrangian (Dolag et 
al. 2005, Vazza et al. 2006; Valdarnini 2010) and Eulerian re- 
simulations of galaxy clusters extracted from large cosmolog- 
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ical volumes found that a sizable amount of pressure support 
(i.e. ~ 10-30 percent of the total pressure inside 0.5/? v „ Q) 
in the ICM can be due to chaotic motions, provided that the 
kinematic viscosity in the innermost cluster region is negligible 
for scales > 10 kpc (Nagai et al. 2007; Iapichino & Niemeyer 
2008; Lau et al. 2009; Vazza et al. 2009; Paul et 2010; Burns, 
Skillman&O'Shea 2010). 

The Adaptive Mesh Refinement technique (AMR) is an optimal 
approach to study the fluid-dynamics of the evolving ICM with 
high spatial resolution, and within a fully cosmological frame- 
work. This technique represents an efficient way to overcome 
the problem of having a too coarse spatial resolution in the cen- 
tral region of collapsed objects, typical of fixed mesh Eulerian 
simulations (e.g. Berger & Colella 1989). The proper trigger- 
ing of the mesh-refinement criteria in AMR simulations allows 
for reaching a peak resolutions comparable to the Lagrangian 
(Smoothed Particles Hydrodynamics, SPH) simulations in the 
innermost cluster regions (e.g., ~ 10-20 kpc), while preserv- 
ing the shock-capturing nature of the algorithm (e.g. Teyssier 
2002). 

Recent works employing cosmological simulations have 
shown that the use of mesh refinement criteria anchored to 
the properties of the 3-D velocity field in the ICM enhances 
the numerical description of chaotic motions (e.g. Iapichino & 
Niemeyer 2008; Vazza et al. 2009; Vazza, Gheller & Brunetti 
2010). Furthermore, the adoption of sub-grid modeling to treat 
the dynamical role played by unresolved gas motions helps to 
model this issue more properly, although the importance of 
these motions in the context of galaxy clusters is not totally 
clear (e.g. Scannapieco & Briiggen 2008; Maier et al. 2009). 

In Vazza et al. (2010, hereafter paper I) we presented a sam- 
ple of 20 massive galaxy clusters re-simulated with the code 
ENZO (Norman et al.2007) and employgin an AMR strategy 
designed to increase the spatial resolution in shocks and turbu- 
lent features in the ICM up to distances of ~ 2 - 3 virial radii 
from cluster centers^. The general properties of the thermal gas 
distribution and of shock waves in the sample were discussed 
in Paper I, in relation to the dynamical history of each cluster 
object. 

In the present paper, we focus on the characterization of the 
chaotic motions of the ICM in these simulated clusters. 
Compared to previous works on the same issue, our approach 
combines a large statistics (which allows us to derive viable 
correlations between turbulent and dynamical features) with 
high spatial resolution (which enables to follow the details of 
the gas dynamics and their spectral properties achieving a good 
scale separation). 

The paper is organized as follows: in Secf2]we present the 
re-simulation technique adopted to produce the sample of clus- 
ters, and its general properties. In Secf3]we discuss the charac- 
terization of turbulent gas motions in the sample, focusing on 

1 In this paper, we adopt the customary definition of the virial ra- 
dius, Ryfr, as the radius enclosing a mean total (gas+DM) overdensity 
of ~ 109 times the critical density of the Universe (e.g. Eke 1998). 

2 A public archive of this data-sample is accessible at 
http://data. cineca. it 
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Fig. 2. Slice showing the curl of velocity for the same region of 
Figffl 



the radial distribution of turbulent motions within clusters, on 
the spectral features of the ICM velocity field across the sam- 
ple and on the scaling relations found between turbulent energy 
and the virial parameters of the host clusters. We explore in de- 
tail the possible connection between the statistical features of 
large-scale turbulent motions in our cluster sample and the ob- 
served statistics of radio halos emission in Sec l3.6l Finally, we 
summarize our conclusions in Sec|4] 

2. Clusters simulations. 

Computations presented in this work were performed us- 
ing the ENZO 1.5 code developed by the Laboratory for 
Computational Astrophysics at the University of California in 
San Diego (http://lca.ucsd.edu). 

The details for the numerical setup are described in Paper 
I. In summary, cosmological initial conditions were produced 
with nested grid/DM particle distributions of increasing mass 
resolution to achieve a high DM mass resolution (wdm = 6.76 • 
10 8 M o ) in the region of formation of each cluster, of volume 
~ (5 - 6R V i r ) 3 , where R V j r is the final virial radius for each 
clusters. 

The clusters were extracted from several boxes sampling a 
total cubic cosmic volume of Lb ox « 440 Mpc hr l . 

We assumed a "concordance" ACDM cosmology with 
Q = 1-0, Qbm = 0.0441, Q DM = 0.2139, Q A = 0.742, Hubble 
parameter h = 0.72 and a normalization for the primordial den- 
sity power spectrum <x 8 = 0.8. 

Our runs neglect radiative cooling, star formation and 
feedback from AGNs, while re-ionization is treated with a 
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Fig. 3. Maps of the absolute value of the chaotic velocity field for three choices of the coherence scale for V\ Umax = 300kpc, 
= 600kpc and = 1200kpc). The side of the slices is 4 Mpc hr l and the depth along the line of sight is 25 kpc h . 



simplified Haardt & Madau (1996) re-ionization model (see 
Appendix of Paper I for details). 

The physical modeling of re-ionization in our runs is 
mainly motivated by the requirement of reproducing realistic 
shock waves, since re-ionization is expected to play a sizable 
role in setting the strength of outer accretion shocks (e.g. Vazza, 
Brunetti & Gheller 2009). 

The clusters in the sample have total masses larger than 
6 • 10 14 M o , 12 of them having a total mass above > 1O 15 M0. 

Grouping the clusters according to their dynamical state is 
useful to detect statistical properties which vary inside the clus- 
ter sample. In Paper I, we outlined our procedure to distinguish 
between "relaxing", "merging" and "post merger" systems of 
our sample, focusing on the total matter accretion history of 
each halo. According to our definition, this sample contains 10 
post-merger systems (i.e. clusters with a merger with a mass 
ratio larger than 1/3 for z < 1), 6 merging clusters at z — 
and 4 relaxing clusters (i.e. systems without evidence of past 
or ongoing major merger for z < 1). 

The refinement strategy inside the w (5 - 6R v , r ) 3 volume 
(hereafter the 'AMR region") was tuned to increase the grid 
resolution based on gas/DM over-density, 5p gas pM and/or ve- 
locity jumps. In this scheme, a normalized 1-D velocity jump 
across 3 close cells in the scan direction (at a given refinement 
level) is computed as 5 W = |Av|/|v m ,„|, where |v m ;„| is the min- 
imum velocity among the 3 cells; the threshold values used 
to trigger the mesh refinement were Sp gas ,DM/Pgas,DM > 3 for 
gas/DM over-densities and 6 V > 3 for the 1-D velocity jump. 
For tests about the numerical convergence of this method and 
for a comparison with the more standard refinement based on 
gas/DM over-density, we address the readers to our previous 
works (e.g. Vazza et al.2009; Vazza 2010). 

The adoption of this refinement criteria was necessary 
to avoid the spurious suppression of turbulent eddies mov- 
ing from dense to less dense regions. This is a major 
drawback of Lagrangian approaches (e.g. Smoothed Particles 
Hydrodynamics, e.g. Dolag et al. 2008 for a review) and of 
grid methods for which AMR is triggered by over-density cri- 
teria only (e.g. Bryan & Norman 1998). In general, the AMR 



region around each of our cluster at z = is always sampled 
with at least N ~ 500 3 cells at the highest available resolution. 

When this extra refinement criterion is adopted, turbu- 
lent eddies can evolve with the maximum possible dynamical 
range, and the decay of turbulence is not artificially damped 
by numerical under-sampling. A first order estimate of the 
Reynolds number available to the largest possible eddies con- 
tained within the AMR region yields a value of the order of 
R e „ jv 4/3 ~ 4000 (e.g. Kritsuk, Norman & Padoan 2006). 
However, based on the velocity power spectra measured for the 
3-D velocity field in our clusters (see Sec l3.3l ), the maximum 
Reynolds number which achieved in our simulated ICM is the 
range of R e ~ N Ali ~ 500 - 1000. In any case, a rather large 
hierarchy of chaotic "eddies" of decreasing size can develop 
within the AMR region and evolve towards the scale of the nu- 
merical dissipation, which is a few times the minimum cell size 

E 

To illustrate the vast amount of spatial details provided by 
the method employed here, in FigQ]we present an example of 
the pattern of temperature fluctuations in a multiple mergers 
event in one of our forming galaxy clusters (El), at the time it 
assembled the bulk of its total mass (z ~ 0.5 - 0.6). The side of 
the slice is « 8.8 Mpc hr x and the width is 25 kpc hr x . The 
rise of fluid instabilities and vorticity at oblique shocks is evi- 
dent behind the wake of accretion shocks at the cluster periph- 
ery, thanks to the sampling at the peak resolution (25 kpc hr l ~) 
even at 4 - 6 Mpc hr x from the cluster center. A patchy multi- 
temperature ICM is observed at this stage of the merger, and a 
multi-scale distribution of eddies is generated in the core clus- 
ter region by the large ~ 1000-2000 km s~ x relative velocities 
of the colliding clumps. 

The small-scale chaotic components of the ICM flow are 
well highlighted by detecting the rotational part of the 3-D 
velocity field (e.g. Ryu et al.2008; Zhu, Feng & Fang 2010; 
Paul et al.2010). Figure [2] shows the map of |V x v\ for the 
same region of FigQ] We compute the vorticity at the scale 



3 For a discussion of the effect of PPM artificial dissipation at the 
smallest scales in the simulation see however the discussion in Secl3.3l 
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Fig. 4. Profiles of the turbulent to total energy ratio for Imax = 
300kpc (top panel) and for Imax = 2n/kMAX (lower panel), at 
z — 0. The different colors refer to clusters with a different 
dynamical state, as explained in Sec|2] 



of the cells with a simple first-order finite difference method, 
in which the difference is computed across a baseline scale of 
Icuri = 50 kpc hT l in the three directions across each cell. We 
find that accretion and merger shocks inject a vorticity of the 
order of ~ 1000- 1500 km s~ l on the scale of l cu ,-i at the location 
of the outer accretion shocks, while shears along the stripped 
infalling sub-clumps (e.g. the cold stream of gas entering on 
the lower right corner in Fig[T]|2| inject small scale vorticity at 
a lower rate, ~ 300 - 600 km s' 1 , but over larger volumes. 



3. Results 

3.1. The decomposition of "bulk" and "turbulent" 
motions. 

In order to characterize turbulent velocity fields in the complex 
environment of galaxy clusters, it is necessary to extract the ve- 
locity fluctuations from the complex 3-D distribution of veloci- 
ties. This is not a trivial task and a number of different strategies 
have been adopted in the recent past to characterize the "bulk" 
and "turbulent" component of the velocity field. This has been 
done for instance by taking spherical shell averages (e.g. Bryan 
& Norman 1999; Iapichino & Niemeyer 2008; Lau et al.2009; 
Burns et al.2010), by mapping in 3-D the velocity field on a 
resolution coarser than the maximum resolution of the simula- 
tion (e.g. Dolag et al.2005; Vazza et al.2006; Vazza et al.2009; 
Valdarnini 2010), or by using sub-grid modeling (e.g. Maier et 
al.2009). 

In this work we adopt two different methods to filter the 
large-scale component of the velocity field of the ICM, in order 
to account for the uncertainty inherent in the particular repre- 
sentation of the turbulent ICM. 

In the first method (hereafter the "A;max" method) we cal- 
culate the power spectrum of the 3-D velocity field of the ICM 
and identified the spatial frequency containing most of the ki- 
netic energy in the flow. According to the simplest view of tur- 
bulent fluids, this maximum scale &max marks the beginning 
of the energy cascade of turbulent eddies down to the dissipa- 
tive scale available to the simulation (which is of the order of 
a few cells at the highest available resolution) and would also 
represent the scale of the maximum Reynolds number of the 
flow. 

Following Vazza et al.(2009) and Vazza, Gheller & Brunetti 
(2010) we measured the power spectrum of the velocity field of 
the simulated ICM, E(k), defined as: 

E(k) = l -m)\\ (i) 

where v{k) is the Fourier transform of the the 3-D velocity 
field (v = [v x , v y , v z ]) defined as : 

Hk) = -i- f v(x)e- Mhx d 3 x. (2) 

The 3-D velocity is measured in the center of (total) mass 
reference frame, and the power spectrum is calculated with a 
standard FFT algorithm on the data at the highest available res- 
olution, with the addition of a zero-padding technique to deal 
with the non-periodicity of the considered volume and of a 
Gaussian apodization function to avoid the generation of spuri- 
ous frequencies at the box edges (see Vazza, Gheller & Brunetti 
2010 and Valdarnini 2010 for a discussion). 

The power spectrum is measured within a cubic region with 
side ~ 4R V j r centered on each cluster, and the scale k^AX is 
set as the scale of the maximum value of k ■ E(k). Then for 
every cluster we filter out the velocity component ~v{k) with 
k < kyiAX, and consider "turbulent" only the inverse transform 
in the real space of the filtered "v{k > kMAx)- For the clusters of 
the sample ky^x is found to correspond to a maximum spatial 
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scale of lyiAx ~ 0.75 - 1.52? V (V ■ Since laminar bulk motions of 
coherence scales smaller than Imax are not filtered out from this 
approach, we expect that this method provides a higher limit to 
the content of the kinetic turbulent energy in the ICM. 

In the second method, following previous works on the 
same topic (e.g. Dolag et al. 2005; Vazza et al. 2006; Vazza et 
al. 2009; Vazza, Gheller & Brunetti 2010) we mapped the 3-D 
velocity field using the fixed spatial scale of Imax - 300kpc as 
a filter to highlight the "turbulent" small scale features of the 
ICM. The filtering procedure works is applied directly in the 
real space, by interpolating the original 3-D velocity field with 
a Triangular Shape Cloud kernel (e.g. Hockney & Eastwood 
1981) of width Imax = 300 kpc to map the local mean field, V\, 
and to detect turbulent fluctuations on small scales. This is mo- 
tivated by early studies based on SPH (Dolag et al.2005) which 
suggested that the typical size of gas/DM clumps crossing the 
virial volume of galaxy clusters with M > \Q ia MqIIi at z ~ is 
< 200 -300 kpc, and that therefore this filtering scale is suitable 
to detect velocity fluctuations on scales smaller than the scale 
of typical laminar infall motion driven by the accreted satel- 
lites. However, in the present work we consider large clusters 
(by a factor ~ 10 in total mass and by a factor ~ 3 in /? Vi >) com- 
pared to previous works. Therefore, a larger physical scale for 
the injection of turbulent motions by accretion of sub-clumps 
may by expected and therefore this may generally yield a lower 
limit on the estimate of turbulent energy within our simulated 
clusters. 

As an example of the effects of adopting different scales to 
estimate the turbulent energy budget, in Fig[3]we show the pro- 
gressive change introduced by the filtering length used to com- 
pute V\ on a cluster simulation. The change of the chaotic field 
with filtering length I is more evident in the outermost regions, 
which are characterized also by bulk motions with coherence 
scales of the order of ~ 300 - 500 kpc (see the evolution of the 
clumps located at the right of the cluster center in FigO. The 
pattern of motions in the central cluster regions are found to be 
tangled at scales much smaller than a Mpc, therefore adopting 
300 kpc as a filtering length (or more) does not produce signif- 
icant differences in the reconstructed turbulent velocity field. 

In the following, we will show that the two filtering tech- 
niques produce statistically consistent results when applied to 
the whole cluster sample, and therefore that the statistical fea- 
tures associated with the turbulent ICM are rather independent 
of the particular filtering method adopted here. 

3.2. Radial distribution of the turbulent energy in the 
ICM. 

The ratio between the turbulent energy of the ICM, E^rb — 
pv^/2 (where v, is the "turbulent" velocity field, v t - v - V\) 
and the total thermal energy, £W™ = 3/2kBTp/(p,m p ) (where 
T is the gas temperature, ks is the Boltzmann's constant, p. is 
the mean molecular mass and m p is the proton mass), is a sim- 
ple and important proxy of the importance played by turbulent 
motions in clusters dynamics. 

Figure [4] shows the behavior of the cumulative ratio 
Euirb I Ethenn inside a given radius for all the clusters in the sam- 
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Fig. 5. Average profiles of the cumulative turbulent to total en- 
ergy ratio at z — 0, for 9 random displacements from the pu- 
tative center of mass of the cluster (solid black profiles). The 
top curves are for a post-merger cluster, the bottom curves are 
for a relaxed cluster. The different colors refer to different dis- 
placements from the center (the corresponding absolute values 
are reported in the panel). 



pie at z — 0; for clarity, we show with different colors the three 
dynamical classes discussed in Secf2] and normalize all radii 
to the R V ir of each cluster. All systems show a profile that in- 
creases with radius, with the smallest values of E Ulr hlE t i ierm (in 
the range < 10-20 per cent for most of clusters) inside 0. 1^ V !>- 
This ratio reaches larger values at 7? v ,>, E tur b/E t i ierm ~ 0.5 - 0.6. 
In the case of relaxed systems the turbulent energy stored in- 
side the core region of the clusters can be as small as ~ 5 
per cent, while in merging or post-merger systems it can be 
as large as ~ 20 - 30 per cent. Inside 0.H? V i>> two systems have 
Eturb/Etherm > 0.3 with both methods, while inside 0.5,R V! > 4 
(7) clusters have Eu, r b / £tfte™ > 0.3 if the filtering is done using 
Imax = 300kpc (Imax = In/kMAx)- 

Some caution should be taken in analyzing the radial dis- 
tribution of turbulence of Fig|4]in the innermost cluster regions 
(r < O.lRyir), since the characterization of the center of mass 
of some systems may be subject to uncertainties, particularly 
in the case of major merger events. In these cases, asymmetries 
in the matter distribution and large (~ 100 kpc) displacements 
between the centers of gas matter and DM can cause uncer- 
tainties in the characterization of the cluster centers. In this re- 
spect, to test the stability of our results, we measure the radial 
distribution of E tur b/E t i ierm (Imax = 300 kpc) assuming 9 dif- 
ferent random centers extracted inside a sphere of w 250 kpc 
around the peak of total mass, for a post-merger cluster and 
for a relaxed one (Figj5j. In both system, the differences in the 
estimated turbulent energy ratio inside r < 0.1/? vlr can be as 
large as a factor ~ 2 - 3; however even in the case of the per- 
turbed system a good convergence in the estimated E tur b/ 'Exh erm 
is achieved for r ~ OAR v( >, even for extreme values of the as- 
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Fig. 6. Average profiles of the turbulent to total energy ratio 
at z = 0, for the 3 dynamical classes of clusters and adopting 
Imax = 300kpc for the filtering of the local velocity field (top) 
or Imax = 2nlk MAX (bottom). 

sumed displacement from the "real" cluster centers. With this 
caveat in mind, in what follows we will mainly focus on the in- 
tegrated values of E Ulr i,IE^ wrm on radii larger than r > 0.1/? V! >, 
observing that in general the differences between the dynami- 
cal classes of clusters are larger than the uncertainty associated 
to the exact location of the cluster center. 

Overall, when compared outside of the innermost re- 
gion, the profiles of turbulent energies in Smoothed Particles 
Hydrodynamics simulations of galaxy clusters (e.g. Dolag et 
al. 2005; Valdarnini 2010) and in AMR grid simulations as 



here (see also Iapichino & Niemeyer 2008; Vazza et al. 2009) 
present a good similarity. However, small sistematic differ- 
ences are present, also due to the different ways in which gas 
matter from infalling clumps is stripped in the two methods 
(e.g. Agertz et al. 2007) and due to the different stratification 
of gas entropy within the inner cluster atmosphere, which in 
turn affects the stability to convective motions in the ICM (e.g. 
Wadsley et al.2008; Mitchell et al.2008; Springel 2010; Vazza 
2010). 

In order to describe more general features of the three 
dynamical classes of clusters, we computed the average pro- 
files of all the cluster within each subsample, as shown in 
Fig|6] Merging, post-merger and relaxing systems define a "se- 
quence", with the latter being less turbulent within the virial 
volume. The larger values of £t«r2> A^tAe™ in merging systems 
suggest that the generation of turbulence starts before the close 
encounter of cluster cores. In Fig|7] we show the evolution 
of gas temperature and total gas velocity for a slice of depth 
25 kpc h~ l through the merging axis of the major merger clus- 
ter E26, for the epochs of z = 0.79, 0.59 and 0.28 (the closest 
cores encounters happens at z ~ 0.3). The actual major merger 
is anticipated by a number of minor mergers of gas/DM mat- 
ter flowing along the massive big filament which connects the 
two clusters, triggering shocks and chaotic motions in the clus- 
ter outskirts. The virialization of the outer ICM is much less 
efficient than in the center, and the temperature is lower com- 
pared to the innermost ICM, therefore it is more easy to drive 
transonic turbulent motions which eventually sink towards the 
center of the post-merger clusters. 

We analyzed the merger sequence of this system by com- 
puting the average values of gas temperature, gas turbulent 
velocity and E Xur blE X h er m in a cylinder of radius « 300 kpc 
along the axis of merger (FigJSJl. Even if at the beginning 
of the merger the turbulent velocity can be fairly small (~ 
50 - 100 km s~'), the energy ratio reaches large values 
(Eturb/E therm ~ 0.4 at z ~ 0.8) due to the presence of cold 
gas falling from the compenetrating outskirts of the two clus- 
ters. During the later stages of the collision, the heating of 
merger shock waves increases the temperature of the ICM and 
decreases the ratio E Ulr b I E t i wrm , even if the turbulent veloc- 
ity field experiences an overall increase (~ 200 km s~ l for 
z < 0.3). This behavior highlights the important fact that the 
ratio between turbulent energy and thermal energy does not al- 
ways mirror a change in the absolute turbulent energy/velocity, 
especially for those regions (or stages) in which the ICM is in 
a highly un-virialized state. 

A qualitatively similar complex picture of galaxy clus- 
ters in a pre-merger phase has also been recently suggested 
by a few X-ray/optical observations (e.g. Murgia et al.2010; 
Maurogordato et al.2010). 

3.3. Spectral properties of the gas velocity fields. 



In Sec l3.ll we introduced the power spectrum of the 3-D ve- 
locity field, E(k), as useful tool to describe the spectral energy 
distribution of the ICM in our simulations. 
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Fig. 7. Time evolution for a slice cut through the center of the major merger cluster E26 (from left to right, images are taken at 
Z - 0.79, z = 0.59 and z = 0.28). The bottom panels show the gas temperature, the top ones show the velocity module. The side 
of each slice is of 8 Mpc h' 1 . 




Fig. 8. Average gas temperature (left), average turbulent velocity (center) and average energy ratio E tllr b / E t i, erm (right) along the 
axis of the major merger of cluster E26, for a cylinder of radius ~ 300kpc. The different lines and colors are referred to different 
redshits. 
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Fig. 9. Average power spectra of the 3-D velocity field for the 
different classes of galaxy clusters in our sample, at z = 0. 



Figure [9] shows the average 3-D power spectra for all clus- 
ters of the sample at z = 0, after averaging each of the 3 
dynamical classes. The spatial frequency k for every cluster 
has been referred to the one corresponding to the virial radius, 
&o ~ 2x/Rvin while E{k) has been normalized to the total ther- 
mal energy inside /? V! >. The power spectra show a well-defined 
power law energy distribution for nearly two orders of magni- 
tude in spatial scale, with a slope of the order of a < 5/3 - 2 
(where E(k) ~ lc~ a ). The maximum of E(k) falls is in the range 
of ~ 1 - 2R V j r and the drop of the power spectra at larger scales 
is clearly detected in our runs thanks to the large volume con- 
sidered. This result basically confirms and extend to larger clus- 
ter masses the results obtained in Vazza et al.(2009) and Vazza, 
Gheller & Brunetti (2010) and shows that the average spec- 
tral distribution of the 3-D velocity field in evolving galaxy 
clusters is self-similar across a wide range of virial masses. At 
spatial scales just below the maximum correlation scale, the 
slope of the power spectra gets significantly steeper compared 
to a standard a — 5/3 slope of Kolmogorov turbulence, while 
it becomes similar to the Kolmogorov behavior at intermedi- 
ate spatial scales. This may be explained as the effect of large 
scale bulk inflows through the cluster virial radius Q, whose ki- 
netic energy adds to the kinetic energy carried by large scale 
turbulent motions. A qualitatively similar picture has been re- 
cently provided also by ENZO-MHD re-simulations of cluster 
formation (Xu et al. 2009, 2010). It should also be noticed that 
the spatial scales smaller than < 8A (where A is the cell size) 
are increasingly affected by the numerical dissipation of the 
Parabolic Piecewise Method employed by ENZO (e.g. Porter & 



4 The motions are likely correlated on large scales due to geometri- 
cal reasons, see Vazza, Gheller & Brunetti (2010). 



Woodward 1994; Kitsionas et al.2009). However, the flattening 
of the spectral slope seen for k/ko > 20 is not dramatic. 

The comparison of the three sub-samples shows that on av- 
erage merging clusters are characterized by significantly larger 
coherence scales. This follows from the fact that, according to 
our definition (Paper I), we identify as merging systems those 
where a close companion is found, in an early infalling phase. 
The larger correlation at ~ 1 - 2/? v „. scales follows from the 
presence of an early stage of interaction between the ICM of 
colliding massive clusters (see also Fig|7]l, whose virial vol- 
umes have just begun to cross each other at the moment of ob- 
servation. 

A complementary tool to measure the spectral features of 
the turbulent ICM is the structure function of the 3-D velocity 
field. As in the case of the Inter Stellar Medium, the structure 
functions (of various order) of the velocity field provide a way 
to compare simulations with the theoretical expectations from 
the basic theory of turbulent flows (e.g. Kritsuk et al.2007). 
In Vazza et al.(2009) we have shown that also in the case of 
the simulated ICM the information provided by the transversal 
and longitudinal structure functions give results consistent with 
the 3-D power spectra, showing maximum coherence scales of 
the ICM flow at ~ /? VJ >, and an overall good consistency with 
standard expectations from the Kolmogorov scaling $3(0 ~ / 
(e.g. Kritsuk et al.2007). 

Here we consider the 3rd order structure function of the 
absolute value of the velocity, |v|: 

5 3 (/)=<||v(r + /)|-|v(r)|| J >. (3) 

For every cluster S 3 was reconstructed by extracting ~ 10 6 
random pairs of cells withing the AMR region and by measur- 
ing the average (volume-weighted) structure function in every 
radial bin. Figure[10]shows the average results for the 3 dynam- 
ical classes of clusters. In this case the distance was rescaled to 
that of the virial radius of every cluster, while the values of S 3 
have been normalized to (where c s is the volume averaged 
sound speed inside each cluster virial radius). The maximum 
of the structure functions is found at the scale of 1 - 2R v ; r , con- 
sistently with the outer scales provided by the power spectra 
analysis. A behavior broadly consistent with a linear scaling, 
S 3(1) ~ l ' 5 ~' is found for all dynamical classes of clusters, 
within the range of 0.05 < r/R vir < 1.5. 

3.4. Scaling laws for the turbulent energy budget in 
clusters. 

The injection of chaotic motions at mergers/accretions is a 
mechanism mainly driven by the gravity of galaxy clusters, and 
therefore overall it should scale with virial cluster parameters. 
One way to model the turbulence injection in the ICM is to 
assume that during mergers a fraction of the PdV work done 
by the sub-clusters infalling onto the main cluster goes into the 
excitation of turbulent motions (e.g. Cassano & Brunetti 2005). 
In this case turbulence is injected in the cluster volume swept 
by the accreted sub-clusters, which is unbound by the effect 
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Fig. 10. 3rd order structure functions for |v(r + /) - v(r)| for 
the dynamical classes of the sample. The spatial scale has been 
normalized for the virial radius of every cluster, while S3 has 
been normalized by (where c s is the volume averaged sound 
speed within R v t r ). 
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Fig. 11. Scaling between the turbulent energy and the thermal 
energy inside a given radius, for three different radii: < /? VI > 
(squares), < 0.5R W i r (stars) and < 0.25/? v „. (crosses). The dotted 
grey line shows the E lur b cc E x i wrm scaling, while the dashed 
gray line shows the E Ulr b ~ 0.27 E t i wrm scaling reported in Vazza 
et al.(2006) for GADGET2 runs. The meaning of colors for the 
data points is the same as in Fig]4] 



of the ram pressure stripping. Since the infalling sub-clusters 
are driven by the gravitational potential, the velocity of the in- 
fall should be ~ V2 times the sound speed of the main clus- 
ter. Consequently, the energy density of the turbulence injected 
during the cluster-crossing should be proportional to the ther- 
mal energy density of the main cluster. Also the fraction of the 
volume of the main cluster in which turbulence is injected (the 
volume swept by the infalling sub-clusters) depends only on 
the mass ratio of the two merging clusters, provided that the 
distribution of the accreted mass-fraction does not strongly de- 
pend on the cluster mass (Lacey & Cole 1993). The combina- 
tion of these points implies that the energy of turbulence should 

3/2 

scale with the cluster thermal energy, E tur b ~ A ■ E t h ~ M 
(where A < 1, e.g. Cassano & Brunetti 2005). In agreement 
with these expectations, Vazza et al.(2006) derived best fit re- 
lations for the scalings between turbulent energy, total virial 
mass and thermal energy in GADGET2 cosmological simula- 
tions, finding E Ulrb oc Mhf and E turb ~ l/3E th . 

In FigQT|we show the scaling between E tur b (measured ac- 
cording to the "kyiAx" filtering) and Exkerm for three different 
radii of integration: r < 0.25/? v »-, r < 0.5/? V! > and r < R W i r . In 
all cases the scaling closely follows E tur b oc Eii wrm ; the scat- 
ter however increases when smaller volumes are considered to 
compute the energies. This supports the idea that the turbulent 
energy measured in our simulations is a fully gravitationally- 
driven process, for which the total mass is the driving param- 
eters. However variations of E tur b/E ( i ierm are found at fixed 
masses and in relation with the dynamical state of the cluster 
(see also Sec l3.5b . 

The turbulent pressure support in the innermost regions of 
galaxy clusters is an important issue because all estimates of 
cluster masses derived by X-ray observables are affected by the 
presence of non-thermal pressure support at some level (Rasia 
et al. 2004; Lau et al. 2008; Piffaretti & Valdarnini 2008; Burns, 
Skillman & O'Shea 2010). Resolved spectra of the chaotic ve- 
locity field are beyond the capabilities of existing X-ray facil- 
ities, at least for the amount of broadening seen by numer- 
ical simulations (e.g. Inogamov & Sunyaev 2003; Dolag et 
al. 2005; Bruggen, Hoeft & Ruszkowski 2005; Vazza, Gheller 
& Brunetti 2010). However, a number of upper limits to the 
amount of turbulent motions in the central region of galaxy 
clusters (e.g. within the cluster cores) are presently available 
(e.g. Churazov et al. 2008; Sanders et al. 2010a). The compari- 
son between these limits and our simulations provides a sanity 
check of our results. 

Very recently, Sanders, Fabian & Smith (2010) published 
upper-limits to the turbulent velocity support available in the 
central region of 62 nearby galaxy clusters, galaxy groups and 
elliptical galaxies. These authors find at least 15 cases with 
less than ~ 20 per cent of the thermal energy density in tur- 
bulence, and weak evidence for turbulent velocities, of the or- 
der of ~ 500 km s~\ in the cool-core cluster RXJ 1347.5-1 145 
that is undergoing a minor merger. It should be stressed that the 
results are obtained after modeling and removing the contribu- 
tion to the broadening of the spatial extent of the objects, that 
is a difficult task (Sanders et al. 2010b and discussion therein). 
Nowadays the most "direct" upper limit is likely that obtained 
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for A1835, which yields E tur i,/E t i, erm < 0.13 on a scale of 
~ 30 kpc (Sanders et al. 2010a). 

Although upper limits to turbulence should be treated with 
caution, in Fig. [12] we plot the upper limits derived from 
Sanders et al. (2010b) for the velocity dispersion versus ICM 
temperature (as vertical arrows) of the 28 galaxy clusters in 
their sample, and the results from our simulations (as open 
squares). The small top squares show the average value of the 
turbulent velocity field ("A;max" method), while the connected 
thick bottom squares show the turbulent energy associated with 
motion on ss 30 kpc scales (which roughly corresponds to the 
projected volume available to the observations of Sanders et 
al. 2010b). Since the turbulent energy spectrum at the smallest 
spatial scales in our simulations may be affected by numeri- 
cal effect in the PPM code at the smallest scales (e.g. Porter 
& Woodward 1994 and Sec l3.31 l. the energy of the turbulent 
motions on scales < 30 kpc was derived analytically from 
the measured total power spectrum on larger scales assuming 
Kolomogorov scaling. 

Figure [12] shows that for most of our clusters the mea- 
sured typical velocity dispersion of motions expected at scales 
I < 30 kpc is of the order of <x,, ~ 0.2c, or smaller, well below 
the upper limits obtained by XMM-Newton. It is interesting to 
notice that, even in the case of merging and post-merger sys- 
tems, where usually very large (I ~Mpc) chaotic eddies with 
Eturb/Etiwrm ~ 0.2-0.3 can develop, the turbulent energy avail- 
able to the smallest scales is consistent with the upper limits 
of Sanders et al. (2010b). We should however notice that ad- 
ditional turbulent injection mechanisms connected with AGN 
outflows are expected to increase the amount of turbulent mix- 
ing around cool core clusters (e.g. Briiggen et al. 2005; Heinz 
et al. 2010; Dubois et al. 2010) and that therefore our estimates 
here can only be strictly compared with non-cool core clusters. 
Importantly, as the authors noted, the upper limits reported by 
Sanders et al. (2010b) may be subject to large uncertainties due 
to the rather uncertain modeling of the thermal gas distribu- 
tion around the target sources. Thanks to the high spectral res- 
olution, in a few years Astro-H will likely provide substantial 
step in our understanding of turbulent motions in the ICM (e.g. 
Takahashi et al. 2010). 

3.5. Time evolution of turbulence. 

The evolution of the total mass for each cluster of our sample 
has been presented in Paper I (Fig. 4). These results have shown 
that the bulk of the total (gas+DM) mass of half of our clusters 
is built in major merger events at z < 1 . This implies that only a 
few dynamical times have past since the major merger forming 
our post-merger clusters. Therefore, correlations between tur- 
bulent energy and the epoch of the last violent merger should 
be expected. 

In order to monitor time-dependent behaviors of turbulence 
in our clusters sample, we measured the turbulent energy (with 
the kuAX method) and the thermal energy also for the additional 
redshifts of z = 0.6 (look back time of « 5.6 Gyr) and z = 0.3 
(~ 3.3 Gyr). Figure [T3l illustrates the evolution of the profile 
of Exurb I Etherm for the clusters classified as relaxed at z = 0; 
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Fig. 12. Scaling between the average temperature and the mean 
velocity dispersion (small squares) for all clusters of the sim- 
ulated sample (squares in colors) and for sample of the clus- 
ter observed with XMM-Newton by Sanders et al.(2010b). The 
additional dotted lines show the dependence of the ICM sound 
speed with the temperature. In order to compare with the obser- 
vations, we filtered the velocity for the same spatial coherence 
scale of I ss 30 kpc available to Sanders et al. (2010b); the data- 
point derived in this way are shown as connected thick squares. 
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Fig. 13. Radial profiles of the E tur b I E t i ierm ratio for the clusters 
classified as relaxed at z = 0. 



the averaging procedure is the same of Sec 13. 41 These systems 
present more chaotic motions at all radii at increasing redshifts; 
in particular at the reference epoch of z = 0.6 their normalized 
profile is on average rather similar to that of major merger sys- 
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Fig. 14. Scaling between the turbulent energy and the thermal 
energy inside a given radius, as estimated for r < R w / r (top 
squares) and for r < 0.257? v( > (lower squares). The dotted gray 
line shows the E Ulr b = Exj wrm scaling, while the long dashed 
one shows the E tur b = Q.21Eii ienn best fit relation of Vazza et 
al.(2006). The meaning of colors is the same as in Fig0] 
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Fig. 15. Energy ratio Eturfc I Ethermal at z = versus the look- 
back time of the last merger event for the post-merger systems 
of the datasample; the horizontal errorbars show the approxi- 
mate uncertainty in the exact epoch of the central phase of the 
major merger event for each cluster. 

support to the merger-turbulence connection in our simulated 
clusters. 



tems at z — 0, suggesting that the ratio E Ulr b/Eii ien „ is strictly 
dependent on the look-back time since the last major merger. 

In Fig[14|we show the values E tur bl E t h erm integrated inside 
spheres of r < 0.25/? v „ ande r < R v / r forz = 0.3 and z = 0.6 for 
all the clusters in our sample (using the "Icmax" method). Also 
in this case the time evolution is rather different when the 3 dy- 
namical classes are compared. Figure[l4]demonstrates the tran- 
sient feature of turbulence in our simulated clusters: variations 
of Exurb/Efaerm ~ 10 are observed in a few Giga-years in con- 
nection with the dynamical status of clusters. Relaxed systems 
move preferentially "downwards" in the (E t b erm , E tur b) plane 
(i.e. by dissipating their turbulent energy and slightly increas- 
ing their thermal energy with time), post-merger systems often 
present a peak of turbulent injection at intermediate redshifts 
(due to the merger event they all experience in the < z < 1 
range), while merging systems present slightly more irregular 
paths across the same plane. A very similar behavior is found 
also by adopting the filter Imax = 300kpc (see also Vazza et 
al.2006). 

Figure Q3] shows the relation between E tllr bl E t i ienn at z — 
(inside r < R w / r ) and the look-back time from the last ma- 
jor merger, tut, for the post-merger systems (adopting the 
"Imax" method). For these clusters we can identify the epoch 
of the last major merger with an uncertainty of the order of a 
cluster-cluster crossing time (~ 1 — 2 Gyr). This Figure sup- 
ports a simple trend between E Ulr bl Eii ierm and tiM, of the type 
Euirbl Etherm 00 —fitiM, with j6 ~ 2/3, and provides additional 



3.6. On the connection between turbulence and giant 
radio halos 

Radio halos are Mpc synchrotron radio sources detected in 
the central regions of a fraction of galaxy clusters (Feretti & 
Giovannini 2000; Kempner & Sarazin 2001; Van Weeren et 
al. 2009; Ferrari et al. 2008 for a modern review). These emis- 
sions are found preferentially in X-ray luminous galaxy clus- 
ters (e.g. Giovannini, Tordi & Feretti 1999; Venturi et al. 2008; 
Giovannini et al. 2009) and only in systems with a X-ray 
perturbed morphology, showing indications of merger activity 
(Buote 2001). Nowadays clear evidence of a statistical connec- 
tion between cluster mergers and radio halos is provided by 
present data (Cassano et al.2010). In principle, several physical 
mechanisms may contribute to the origin of non thermal com- 
ponents to power such large scale radio emission (e.g. Blasi 
& Colafrancesco 1999; Sarazin 1999; Dolag & Ensslin 2000; 
Brunetti et al. 2001) 

Turbulence in the magnetized ICM may power large scale 
radio emission in galaxy clusters, provided that relativistic 
electrons can couple with the MHD modes excited during 
merger events (e.g. Brunetti et al. 2008 and references therein). 
Brunetti & Lazarian (2007, 2010) developed a comprehensive 
modeling of the properties of MHD turbulence in the ICM and 
of turbulent acceleration of relativistic particles. According to 
their picture compressible turbulence provides the most im- 
portant contribution to the process of particle (turbulent) re- 
acceleration in the ICM. Under the assumption of a collisional 
coupling between the turbulent modes and both the thermal and 




Fig. 16. Cumulative volume distribution functions for the ratio between turbulent energy and thermal energy inside ~ IMpc 3 for 
two clusters in the sample, by adopting the filtering with Imax = 300kpc. 



relativistic particles in the ICM, these authors have shown that 
radio halos can be switched on in massive clusters when com- 
pressible turbulence (generated on ~ 300 kpc scales) accounts 
for > 15-30 per cent of the thermal energy. Although we 
have shown in Sec l3.2l that in merging and post-merger clus- 
ters E tur b/E t herm is fairly large, turbulence should also be rather 
volume filling in order to produce radio halos of ~ Mpc size. 

In Fig[l6]we show the low redshift (z < 0.2) evolution of 
the cumulative volume distribution of E tur b/ 'E therm, inside a ref- 
erence fixed volume of the order of that filled by typical radio 
halo IMpc 3 ) centered in the cluster center of two massive 
systems of our sample, with nearly identical final mass: a re- 
laxed cluster at z = (right) and a major merger cluster (left). 
To be conservative in Fig[l6]we used the Imax = 300 kpc filter, 
which likely provides a lower limit to turbulent motions in the 
ICM. 

The fraction of volume where E Ulr bl E t i ienn is above a fixed 
threshold is larger in the post-merger cluster. At z ~ only ~ 5 
per cent of this central volume has E^b/Etherm > 0.3 in the 
relaxed cluster, compared to a ~ 30-40 per cent in the post- 
merger cluster. Although turbulence within the scale length of 
< 300 kpc is not totally volume filling even in the post-merger 
cluster, the whole total "projected" volume within the central 
1 Mpc 3 of the cluster would be seen as "turbulent" from ev- 
ery line of sight (e.g. Subramanian et al. 2006; Iapichino & 
Niemeyer 2008). Similar results can be found for all the clus- 
ters within the sample, showing a very tight correlation be- 
tween the volume filling factor of turbulence and the dynamical 
state of the host galaxy cluster. 

Our sample is large enough to allow a first statistical esti- 
mate of the frequency of "turbulent" clusters and a first compar- 
ison with the observed frequency of radio halos in real galaxy 
clusters. The cluster-mass distribution of completeness of our 



sample is consistent with the Press & Schecter mass function 
for M > 7 • 1O 14 M0//i; for lower masses, the lack of objects is 
of the order of ~ 30 per cent (FigfTTli. This shows that our sam- 
ple can be safely used for statistical studies provided that our 
clusters are not biased with respect to a particular dynamical 
state (as in our case, see Paper I). 

According to calculations of turbulent acceleration in 
galaxy clusters an energy ratio of E mr b I Etherm « 0.2 • 
y/S/iksTkev) (where T^ e v is the cluster temperature measured 
in keV) may allow the generation of radio halos (Brunetti 
& Lazarian 2007) B In Fig J 1 81 we show the distribution of 
Eturb/E t i ierm for the clusters in our sample, evaluated in the typi- 
cal volume occupied by radio halos. The size of observed radio 
halos is usually in the range of ~ Mpc, that corresponds to 
~ R V ir/3 for most of the clusters of our sample (for complete- 
ness we also report the distribution for r < /? v ,>/10 and R V i r /2 
in the same Figure). The cumulative distributions extracted for 
< R V ir/2 and /? v ;,73 are similar for the Imax = 300kpc fil- 
tering (right panel), while the turbulent energy ratio increases 
in the Imax filtering when larger volumes are considered for 
the integration; this may be explained with the large scale 
bulk motions present outside of the cluster innermost region. 
Regardless of the filtering technique, ~ 1/3 of the massive 
galaxy clusters in the local Universe is found to host turbulent 
motions of the order of E t herm/4 within R V j r /3- This suggests 
that the theoretical conditions for the generation of radio halos 
in our simulated clusters are achieved for ~ 1/3 of our clusters, 
provided that turbulence drives the formation of diffuse radio 
emission. Nowadays extensive radio follow ups of complete 
samples of galaxy clusters indicate that a ~ 1/3 of massive 
L x > 5- lO^erg/s (M > 10 15 Mq) objects host Mpc-scale halos 

5 By considering compressible turbulence generated at ~ 300 kpc 
scale. 
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Fig. 17. Cumulative mass function for the clusters of our sam- 
ple (solid line) compared to the theoretical Press-Schecter mass 
function for the same volume (dashed line). The vertical error 
bars mark the Poissonian noise for every mass bin. 



(e.g. Venturi et al.2007; 08; Cassano et al.2008) and that all ra- 
dio halo clusters are merging systems (e.g. Cassano et al.2010). 
This is consistent with our previous claim. Semi-analytical cal- 
culations based on extended Press-Schechter theory already 
demonstrated that the theoretical occurrence of Mpc-scale ra- 
dio halos (assuming that they originate from turbulent accel- 
eration) is consistent with the occurrence of these sources as 
measured by present surveys (Cassano et al.2008). This is the 
first time in which an overall consistency between the fraction 
of "turbulent" clusters and that of radio-halos clusters in mas- 
sive objects is underlined by means of numerical cosmological 
simulations. 



4. Discussion and conclusions 

In this paper we explore the properties of turbulent motions in a 
statistical sample of massive galaxy clusters simulated at high- 
resolution with the AMR code ENZO (Norman et al.2007). 

We focus on a sample of 20 massive galaxy clusters already 
presented in Paper I and simulated with the tailored AMR tech- 
nique introduced in Vazza et al.(2009). We provide a complete 
view of the onset and of the evolution of gravity-driven turbu- 
lent motions which are injected in massive galaxy clusters via 
mergers. 

However, the physical modeling adopted in our simulations 
does not include the effect of AGN activity, Cosmic Rays and 
magnetic field. In reality, AGN activity may drive small-scale 
turbulent motions around the core region of clusters (e.g. Heinz 
et al.2006; Scannapieco & Bruggen 2009; Dubois et al.2010, 
Morsony et al.2010; Vazza 2010), while the interplay of tan- 
gled magnetic field and Cosmic Rays is expected to inject 
chaotic motions at very small scales, via instabilites in the ICM 
(Parrish & Stone 2008; Quataert 2008; Ruskowsky & Oh 2010; 
McCourt et al.2010; Ruszkowsky et al.2010). 



The magnetic field is also expected to play a role along 
the cluster evolution (e.g. Dolag 2006; Dubois & Teyssier 
2008;Dolag & Stasyszyn 2009; Donnert et al. 2009; Collins et 
al. 2010), even if without affecting the spectral properties out- 
lined in Sec l3.3l for a wide range of scales (e.g. Xu et al.2009; 
Xu et al.2010). 

The implicit assumption to trust in the development of the 
turbulent motions as derived from our simulations is that kine- 
matic viscosity in the ICM is negligible for all scales larger 
than the minimum cell size, 25 kpc h~ l , which implies that 
the effective kinematic viscosity, v, must be smaller than v ~ 
10 29 cm 2 s _1 [§. At smaller scales there are physical reasons to 
believe that the turbulence in the ICM is finally dissipated in a 
collision-less regime, e.g. by accelerating relativistic particles 
(Brunetti & Lazarian 2007). A correct modeling of the dynam- 
ics of turbulent motions at the scales between the numerical 
resolution of the scheme and the physical scale responsible for 
dissipation (which is likely << kpc) can be computed by means 
of sub-grid modeling, by assuming some closure relation for 
the turbulent equations. Valuable attempts have been recently 
done in the framework of galaxies and galaxy clusters simu- 
lations (e.g. Scannapieco & Bruggen 2008; Maier et al.2009; 
Iapichino et al.2010); however the application of these tech- 
niques in astrophysics is still in its beginnings. 

To summarize our results, the analysis of turbulent motions 
in this sample of massive galaxy clusters extends the previous 
results based on the same AMR technique in ENZO and fo- 
cused on clusters with smaller masses (Vazza et al.2009; Vazza, 
Gheller & Brunetti 2010). Given the fairly large number of ob- 
jects in this present sample (20) we were able to study the de- 
pendence of the turbulent features on the dynamical state of 
clusters, following from their matter accretion history across 
cosmic time. Two methods were presented to detect turbu- 
lent motions in the ICM. One is based on a filtering in the 
Fourier space of the component of velocities associated with 
wave numbers larger than the wavenumber of the maximum 
spectral energy, and the other is based on the filtering in the 
real space of the velocity component with coherence scales 
smaller than the fixed length of Imax = 300 kpc. (Sec B.U . 
We have shown that the statistical results obtained in this pa- 
per are largely independent of the particular method adopted. 
Post-merger and merging clusters show large values of turbu- 
lent energy compared to the thermal energy of the ICM, with 
E Ulr hl E t i ienn ~ 0.2 - 0.3 in the innermost cluster regions. On 
the other hand relaxed clusters show much lower values of the 
turbulent ratio, E tur b/E t f, erm ~ 0.05 (Sec l3.2b within the same 
radius. These results are in line with recent studies performed 
by Paul et al. (2010), but extends to clusters with larger masses. 

Owing to the very high dynamical range achieved in our 
simulations (N ~ 550) we were able to study the spectral fea- 
tures of the ICM velocity field with an unprecedented separa- 
tion between the forcing and the dissipation scales, achieving a 



6 This can be roughly estimated considering that at the scale of the 
minimum cell size, Ax = 25 kpc, the typical value of the disper- 
sion in the velocity field is v < 50 km s , as shown in Fig[l2] and 
the Reynolds number is equal to one, and therefore : v ~ 50km s • 

25 kpch- l /R e ~ 10 2< W s-\ 
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Fig. 18. Distribution functions for E^b/Etherm inside three reference radii (R v( >/2, R V i r /3 and i? vi >/10) for the simulated clusters 
at z — 0. The left panel is for the turbulent ICM velocity field estimated using the Imax filtering, while the right panel is for the 
filtering at Imax = 300kpc. The solid lines show the differential distributions (for r < 7? vl >/3 only), while the dashed lines show 
the cumulative distributions. The vertical grey band shows the approximate regime of turbulence required by present calculation 
ofr the turbulent re-acceleration of relativistic electrons (e.g. Brunetti & Lazarian 2007). 



typical Reynolds number of the order of R e ~ 500 - 1000. The 
power spectra of the 3-D velocity fields extend across nearly 
two orders of magnitude, with E(k) ~ k '^~ 2 , and show a typ- 
ical scale for the peak of the energy spectrum at the scales of 
1 - 2R V j r (with the tendency of merging clusters to present the 
largest outer correlation scales). Consistent results were also 
obtained by the measure of the third-order structure function 
of the velocity field across the clusters sample (Sec l3.3l l. As a 
sanity check we compared our results with the available lim- 
its on turbulent motions in the ICM for X-ray observations 
(Sanders et al. 2010b). We show that the available limits on 
turbulent motions on ~ 10-50 kpc scales obtained in the com- 
pact cores of relaxed clusters, E tl , r i,/E t i, erm < 10 per cent, are 
well consistent with the amount of these motions measured in 
our simulated clusters (Sec l3.41 >. We derive scaling laws for the 
integrated turbulent energy in our clusters, and confirm previ- 
ous results based on GADGET2 runs (Dolag et al. 2005; Vazza 
et al. 2006), where turbulence is found to scale with the ther- 
mal energies (Sec l3.4t . The time evolution of the turbulent en- 
ergy within clusters was sampled at different epochs (z = 0, 
Z = 0.3 and z = 0.6), reporting a significant similarity between 
the turbulence found in the past of relaxed clusters and that 
found in merging clusters at z — (Sec l3.51 >. We also found 
an anti-correlation between the look-back time since the last 
major merger of a cluster and the level of turbulence at z = 0. 

Finally we explore the connection between turbulence and 
radio halos in galaxy clusters. Thanks to the unprecedented 
statistics of our sample of massive clusters simulated with 
AMR, we can compare the occurrence of turbulence in our 
clusters with that of observed giant radio halos in nearby mas- 
sive clusters. Current calculations of turbulent acceleration for 



the origin of radio halos suggest that these sources can be gen- 
erated in the turbulent ICM (e.g. Brunetti & Lazarian 2007). 
We found that these conditions are reached only in the case of 
merging systems, where ~ 1 /3 of the cluster volume is in the 
form of turbulent motions for a few Gyr, while only a few per- 
cent of the cluster volume is turbulent in the case of relaxed 
systems (Sec l3.61 l. In particular we find that ~ 1/3 of our clus- 
ters show E tur b/E t f, em , > 0.15 - 0.30 level on Mpc-scale and 
we notice that this fraction is consistent with that of massive 
clusters hosting giant radio halos. These findings are in line 
with previous results based on extended Press-Schechter the- 
ory, that show that the theoretical occurrence of radio halos in 
the turbulent acceleration model is consistent with present ob- 
servations. 



The generation and evolution of turbulent motions in our 
simulated ICM is very complex, and it follows from the multi 
scale energy distribution of the turbulent ICM, as produced by 
the dynamical evolution of galaxy clusters. It is very intriguing 
that our rather simple physical setup (which does not consider 
any non-gravitational mechanism) is sufficient to capture some 
of the most important multi scale features at work in the tur- 
bulent ICM and to produce results in line with existing obser- 
vations. The adoption of an ad-hoc Adaptive Mesh Refinement 
scheme specially tailored is mandatory to capture the full range 
of scales needed to describe the evolution of turbulence in the 
ICM; however in the future more sophisticated physical mod- 
eling including magnetic fields,Cosmic Rays would be needed 
to model also the morphological and spectral features of syn- 
chrotron radio emission in a detailed and time-dependent way. 
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